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ABSTRACT 

In a recent analysis of the 7-year WMAP temperature sky maps, Gurzadyan and Penrose claim 
to find evidence for violent pre-Big Bang activity in the form of concentric low-variance circles at 
high statistical significance. In this paper, we perform an independent search for such concentric low- 
variance circles, employing both statistics and matched filters, and compare the results obtained 
from the 7-year WMAP temperature sky maps with those obtained from ACDM simulations. Our 
main findings are the following: We do reproduce the claimed ring structures observed in the WMAP 
data as presented by Gurzadyan and Penrose, thereby verifying their computational procedures. How- 
ever, the results from our simulations do not agree with those presented by Gurzadyan and Penrose. 
On the contrary we obtain a substantially larger variance in our simulations, to the extent that the 
observed WMAP sky maps are fully consistent with the ACDM model as measured by these statistics. 

Subject headings: cosmic microwave background — cosmology: observations — methods: statistical 



1. INTRODUCTION 

One of the main accomplishments in cosmology during 
the last decade is the establishment of the ACDM in- 
flationary concordance model. According to this model, 
the universe consists of 5% ba ryonic matter 22% dark 
matter and 73% dark energy (jJarosik et al.ll20l"o|) . and 
is filled with random and Gaussian fluctuations. These 
fluctuations were generated during a short period of ex- 
pone ntial expansion called inflation (e.g. Liddle fc LvthI 
120001 and references therein) , during which the universe 
expanded by a factor of ~ 10^^ in ~ 10"'^^ s. With only 
a handful of free parameters, this model is able to suc- 
cessfully fit thousands of observational data points. 

Nevertheless, the ACDM model must at the current 
stage be considered an effective model rather than a fun- 
damental model. First, it relies on several quantities 
that have never been directly observed except through 
their gravitational impact, such as both dark matter and 
dark energy. Second, it postulates the existence of an 
unknown scalar field, the inflaton. It is therefore impor- 
tant to put the inflationary framework to stringent tests, 
probing its range of validity in different ways. One ap- 
proach to do so is to construct alternative cosmological 
theories, making different observational predictions than 
ACDM, and then compare the two models using high- 
precision data. 

One example of such work has been demonstrated by 
the development of the "Conformal Cyclic Cosmology" 
(CCC) model bv IPenrosd ((20081 [20091 [2010l) . In this pic- 
ture, the history of the universe is described in terms 
of a series of "aeons" , each of which is defined as a fi- 
nite period between two Big Bang events. We will not 
consider further details of the CCC model in this pa- 
per, except for one interesting feature: According to 
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IGurzadvan fc Penrose! (j2010l) , the model postulates that 
super-massive black holes may collide in earlier aeons, 
releasing tremendous amounts of energy in the form of 
gravitational waves. These collisions may be observable 
in our current aeon in the form of concentric circles of low 
variance in the cosmic microwave background (CMB). 

Foll owing up on this prediction, IGurzadvan fc Penrose! 
(|2010f) analyzed the 7-year WMAP temperature sky 
maps, searching for concentric circles of low variance. 
And quite surprisingly, they claim to find evidence for 
such rings at the 6cr confidence level, by comparing the 
results obtained from WMAP with ACDM simulations. 
This claim is sufficiently spectacular (involving both pre- 
Big Bang phenomena and super-massive black holes) to 
catch the interest of the general media, with numer- 
ous news stories and live media appearances following 
in the weeks after the release of the paper. Given these 
widespread reactions, the claims of Gurzadyan and Pen- 
rose deserve closer scrutiny through independent analy- 
sis; this paper presents one such independent analysis. 

2. DATA AND SIMULATIONS 

The pre-Big Bang b l ack h ole collisions discussed by 
'Gurzadva n &: Penrose! (|2010l ) observationally manifest 
themselves in the form of one or more rings of low 
variance centered on a central position, p, in the CMB 
anisotropy field, AT(n), observed today. In this paper, 
we therefore consider the best currently available full-sky 
maps of the CMB, namely the 7-year WMAP tempera- 
ture sky maps. These data are provided in the form of 
pixelized HEALPix sky maps with a pixel resolution of 
7'. 

For simplicity we consider only the foreground-reduced 
WMAP W-band (94 GHz) data in this paper, as 
this channel has the highest resolution of all the 
WMAP bands, with a beam of 13' FWHM. (Note that 
IGurzadvan &: Penrose! (^010) additionally analyses the 
WMAP V-band data and the BOOMERanG98 data; 
however, we agree that instrumental effects are irrele- 
vant to this analysis, and therefore consider only the W- 
band in the following.) We apply the WMAP KQ85 
sky cut to this map, removing 22% of the sky, and 
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Fig. 1. — Examples of single standard deviation profiles (thick histograms), computed from WMAP (left) and a simulation (right). The 
thin solid line shows the mean of all profiles, and the shaded regions show the 1 and 2u confidence regi ons. The WMAP profile is centered 
on galactic coordinates (Z,b) = (105.04°, 37°), reproducing Figure 2 of Gurz advan fc Penros"3 120101) . Note that low-variance rings are 
found also in the simulation. 



taking into account bo th the galactic pl ane and high- 
latitude point sources (IGold et all [20T1 . (Note that 
IGurzadvan fc Penrosi ()201ClO adopts a straight |6| < 20° 
sky cut; the difference is not significant.) 

Next, we build an ensemble of 1000 simulations with 
a spect rum given by the be st-fit WMAP 7-year ACDM 
model (jKomatsu et al.|[201Cl[) . These realizations are con- 
volved with the W-band beam and HEALPix pixel win- 
dow, and projected onto a HEALPix A"sido = 512 grid. 
Finally, uncorrelated Gaussian noise with RMS equal to 
fo / V A^obs (p) are added to pixel p, where (Tq is the W- 
band noise RMS per observation and A'obs (p) is the num- 
ber of observations in that pixel. The KQ85 sky cut is 
also applied to the simulations. 

3. METHOD 

The central concept in the following search is the radial 
standard deviation profilcQ, 
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where p denotes a fixed reference pixel on the sky, n is 
a variable pixel index, and 9 — acos(p • n) is the angu- 
lar distance between those two positions; b denotes an 
angular bin, [0f,_i,0f,)- Thus, this function measures the 
standard deviation of the CMB field as a function of ra- 
dial distance from p, and a sharp negative spike in this 
function corresponds to the postulated signature of the 
black hole collisions. To reduce correlations between dif- 
ferent bins, we subtract the mean from each function 
crp(b) separately before further processing. 

We measure ap{b) for all positions over a HEALPix 
A'sidc = 32 grid, for a total of 12 288 points. Each func- 
tion is evaluated from to 20° in 40 bins, corresponding 
to a bin size of 0.5°. If more than 40% of the pixels 
within 20° of a given pixel p are masked out, that pixel 
is removed from further analysis. Finally, we also remove 
any pixels that are immediate neighbors to the galactic 
cut. 

* We use the terms "standard deviation profile" and "variance 
profile" interchangeably in this paper. 



Having computed these functions for both data and 
simulations, the next step is to identify potential low- 
variance ring candidates, and quantify their significance. 
To this end, we introduce two different statistics. First, 
we consider matched filters designed to highlight negative 
spikes in <Tp(b), 



(2) 



Here f{h) defines a discrete filter, and we consider 
three different cases in this paper, namely /i = 
[0.5,-1,0.5], /2 = [0.25,0.25,-0.5,-0.5,0.25,0.25] and 
h = [0.5,-0.125,-0.75,-0.125,0.5]; for \h-h'\ larger 
than the length of the filter, / is zero. These filters cor- 
responds to 1) a negative top-hat filter of width 1; 2) 
a negative top-hat filter of width 2; and 3) a negative 
wedge filter, each sensitive to typical interesting candi- 
dates. For each case, we adopt the maximum of dp{b) as 
our statistic, considering only angular distances larger 
than 2° as the intrinsic estimator variance is very large 
on the smallest scales. We then make both sky maps and 
histograms of the maximum values of (7p(6), and compare 
these between the observed data and the simulations. 

Our second statistic is based on the standard esti- 
mator, which is sensitive to the overall fluctuation level 
of (Tp{b), rather than individual spikes. Specifically, we 
compute 
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for each pixel and data set, where fib is the mean of 
(Tpib) computed from the simulation set, and Cbb' is the 
covariance matrix. Again, we make sky maps of these 
values for the observed data, and compare the observed 
data with the simulations in terms of histograms. 

4. RESULTS 

First, we show in the left panel of Figure [1] 
the standard deviation profile computed from the 
WMAP data centered on Galactic coordinates {l,b) — 
(105.04°, 37°). (The mean is not subtracted in this 
plot.) This is the same profile as shown in Figure 2 of 
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iGurzadvan fc Penrose! (|20lol ). The agreement between 
the two results is excellent, despite the different masks 
used by the two analyses. This validates the routines 
used to compute the radial profiles for both analyses. 

We have also plotted the mean of the variance profiles 
computed from all pixels in Figure [1] together with the 
corresponding 1 and 2cr confidence regions. Here we note 
several interesting features. First, the WMAP example 
profile is consistently low compared to the mean, sug- 
gesting strong correlations between bins. This is typical 
for all profiles for both WMAP and simulations; the local 
variance in a CMB map depends on the scanning strat- 
egy of WMAP, and the entire profile can shift up or down 
depending on whether the corresponding pixel is located 
in the ecliptic plane (high noise) or in the ecliptic poles 
(low noise). Since this effect is of little interest in the 
search for concentric rings, we subtract the mean before 
further analysis to reduce bin-to-bin correlations. 

Second, the simulated distribution is very similar to 
that of WMAP, having consistent mean and standard 
deviation. We have also checked that this distribution is 
consistent from simulation to simulation. 

Third, we see that the lowest point in the WMAP 
example corresponds to a ~ 3.3(t outlier compared to 
the typical WMAP profile. This is the point that was 
claime d to be a 6<7 outlier by IGurzadvan fc Penrose! 
((Mol) . With our simulations, this point appears consid- 
erably less anomalous. Further, it is important to note 
that this particular profile is by no means randomly cho- 
sen. On the contrary, it was picked out precisely due 
to its anomalous behaviour after searching through more 
than 10000 points. The actual statistical significance of 
this particular B.Scr outlier is therefore likely quite low. 

Fourth, as illustrated by the example simulation profile 
in the right panel of Figure [U it is not difficult to find in- 
dividual simulated profiles with "anomalous" behaviour. 
Here one can also see "rings", similar to those observed 
in WMAP, but created by chance alone. 

To address these issues more rigorously, we are forced 
to adopt statistical techniques; visual inspection is not 
sufficient. We therefore employ the search algorithms 
described in Section |3] to identify peculiar candidates 
in both the data and simulations, and consider various 
statistics based on these results to interpret significances 
statistically. 

In Figure[2l we show sky maps of the third matched fil- 
ter, /3, statistic, where each pixel indicates the maximum 
of (Tp(6) over b. Thus, pixels with large values in these 
plots indicate variance profiles with a notable negative 
wedge-like structure. Such maps are shown for both the 
WMAP data and a random simulation. At least visually, 
the two maps appear statistically consistent. 

This statement is quantified more rigorously in Figure 
[3l where we show histograms for each of the four statis- 
tics computed from the WMAP data, and compared to 
the mean properties of the simulated ensemble. Again, 
the WMAP properties appear fully consistent with the 
simulations. 

However, the histogram statistics shown above are sen- 
sitive only to the mean properties of the CMB field. They 
therefore only show that there is no evidence for a large 
number of concentric circles in the WMAP data, not 
that there cannot be a small number of highly significant 
cases. Finally, we therefore also consider the maximum 




Fig. 2. — Example of sky maps obtained by the matched filter 
searches, in this case the wedge filter, /s. The top panel shows the 
map obtained from WMAP, the bottom panel shows the same for 
a random simulation. 
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Fig. 3. — Histogram comparison between the WMAP data (thick 
solid lines) and the ACDM simulations (mean in thin black lines; 
1 and 2(T regions in gray bands), for each of the four statistics 
adopted in this paper. 

value of each of the four statistics, which should be sen- 
sitive to single extreme cases. Figure |4] shows histograms 
of such maximum values, as computed from the sim- 
ulated ensemble, with the corresponding WMAP value 
indicated by a vertical line. Again, we see that WMAP 
appears fully consistent with the ACDM simulations. 

5. CONCLUDING REMARKS 

In this paper we search for concentric low-variance 
rings in the 7-year WMAP temperature sky maps, 
seeking to reproduce t he re sults recently presented by 
IGurzadvan fc Penrose! ()2010f) . While our two analyses do 
agree in terms of specific variance profiles, they clearly 
disagree in terms of statistical interpretation and signifi- 
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Fig. 4. — Comparison of the full-sky maximum value for each 
statistic between the simulations (histograms) and the WMAP 
data (vertical dashed line). 

cance. Specifically, we find a subst antially larger variance 
in our ACDM simulations than iGurzadvan fc Penrose! 
([2010) do in theirs. When taking into account this larger 
variance, and accounting for statistical selection effects, 
the evidence for concentric circles in the WMAP data ap- 
pears minimal. Rather, the WMAP data appears fully 
consistent with our simulations. 

The main difference between the two analyses must lie 
in the construction of the simulations. However, we have 
good reason to believe that our simulations are correct. 
First, we note that in our simulations the mean variance 
profile decreases towards small angular distances. This is 
bound to happen for two reasons: The CMB anisotropics 
constitute a correlated field with a characteristic scale of 
1°, corresponding to the first peak in the CMB spectrum. 



Also, the instrumental beam of WMAP smooths out all 
small scale structure. It is theref ore surprising that this 
effect is not seen in Figure 2 of IGurzadvan fc Penrose! 
(|2010f l. Second, the fact that our simulations agree 
with WMAP provides further confidence; it is difficult 
to imagine an error in the simulation pipeline that would 
cause the simulations to become more similar to the ob- 
served data. Third, as noted by IGurzadvan fc Penrose! 
(|201Cll) , there are both low and high peaks in the reported 
variance profiles. While Gurzadyan and Penrose assign 
no significance to the high variance peaks (stating that 
"the peaks of high variance are of no importance, as these 
can result from numerous irrelevant effects"), they are 
clearly relevant in our interpretation, as they illustrate 
the substantial statistical variance present in these pro- 
files. This large variance is observed both in the real data 
and the simulations. 

Thus, we conclude that there is no evidence for the 
CCC m odel in the current WMAP data. Of course, 
Planck (jTauber et al.! I2O1O0 may provide new light on 
this issue, having higher sensitivity and resolution than 
WMAP, but one should probably not have too high ex- 
pectations in this regard. Even if the CCC model should 
turn out to describe the real universe, it will quite likely 
be difficult to unambiguously identify such concentric cir- 
cles due to the dominant background cosmic variance. 

The computations presented in this paper were car- 
ried out on Titan, a cluster owned and maintained by 
the University of Oslo and NOTUR. Some of the re- 
sults in this paper h ave been derived using the HEALPix 
(jGorski et al.l 120051 ) software and analysis package. We 
acknowledge use of the Legacy Archive for Microwave 
Background Data Analysis (LAMBDA). Support for 
LAMBDA is provided by the NASA Office of Space Sci- 
ence. 
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